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ABSTRACT 



We present the first ground-based detection of thermal emission from an exoplanet in the //-band. Using HAWK-I on the VLT, 
we observed an occultation of WASP- 19b by its G8V-type host star. WASP- 19b is a Jupiter-mass planet with an orbital period of 
only 19 h, and thus, being highly irradiated, is expected to be hot. We measure an //-band occultation depth of 0.259^ ^ %, which 
corresponds to an //-band brightness temperature of T H = 2580 ± 125 K. A cloud-free model of the planet's atmosphere, with no 
redistribution of energy from day-side to night-side, under predicts the planet/star flux density ratio by a factor of two. As the stellar 
parameters, and thus the level of planetary irradiation, are well-constrained by measurement, it is likely that our model of the planet's 
atmosphere is too simple. 
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1. Introduction 

For a planet that is occulted by its host star, we can measure the 
planet's emergent flux and thus determine its brightness temper- 
ature without spatially resolving the system (e.g., Charbonneau 
et al. 2005; Deming et al. 2005). From this, we can determine 
a planet's albedo and the efficiency with which energy is re- 
distributed from the day-side to the night-side of the planet 
(e.g., Barman et al. 2008; Alonso et al. 2009a). By measur- 
ing occultations over a range of wavelengths, we can map the 
planet's spectral energy distribution (SED) and can thus infer 
its chemical composition (e.g., Swain et al. 2009a), investigate 
whether a temperature inversion exists (i.e., a stratosphere; e.g., 
Charbonneau et al. 2008; Knutson et al. 2008), and constrain the 
albedo and the day/night energy redistribution more strongly. 

The Spitzer Space Telescope (Werner et al. 2004) has mea- 
sured planetary thermal emission in the 3.6-24ji/m range for 
15 exoplanets (e.g., Deming et al. 2009; Cowan & Agol 2010; 
and references therein). These observations appeared to iden- 
tify two classes of short-period, giant planets, those with strato- 
spheres and those without, depending on whether there are 
upper-atmosphere absorbers, and hence the planet's effective 
temperature (e.g., Fortney et al. 2008). Retarded cooling pro- 
duced by opaque atmospheres is one possible cause of the 
anomalously large measured radii (compared to the predictions 



* Based on data collected with the VLT/HAWKI instrument at ESO 
Paranal Observatory, Chile (programs 083.C-0377(A)). 
** The photometric time-series used in this work are only available 
in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 



of simple models) of several exoplanets (e.g., Burrows et al. 
2007). 

For some inflated planets, such as WASP- 17b (Anderson et 
al. 2010), tidal heating from the circularisation of an eccentric 
orbit seems a more likely explanation. Measuring orbital eccen- 
tricity is essential to constrain the current and past rates of tidal 
heating, and the timing of occultations is essential for constrain- 
ing the eccentricity. 

The Spitzer measurements probe only the uppermost layers 
of the exoplanet atmospheres studied to date because their SEDs 
peak in the 2-4 pm range. Observations at shorter wavelengths 
are needed to probe the deeper layers of the atmosphere where 
energy redistribution takes place (Burrows et al. 2007), and to 
provide greater leverage when trying to estimate the slope of 
the temperature structure. Giant planets with orbital periods of a 
day and shorter are starting to be discovered (Hebb et al. 2009; 
Hellier et al. 2009; Hebb et al. 2010). With SED peaks around 
1 pm, measurements bluewards of Spitzer wavelengths are even 
more important for these planets to estimate the bolometric lu- 
minosity and thus measure the temperature, albedo, and energy 
redistribution efficiency. 

Spitzer 's wavelength coverage is invaluably extended by 
ground-based facilities, which have measured occultations in the 
near-infrared (near-IR) detect: OGLE-TR-56b (z'-band, Sing & 
Lopez-Morales 2009), TrES-3b (/T-band, de Mooij & Snellen 
2009), and CoRoT-lb (2.09 pm, Gillon et al. 2009a; and K s - 
band, Rogers et al. 2009). The Hubble Space Telescope has also 
measured occultations in the near-IR for the two exoplanets with 
the brightest host stars (HD 189733b, Swain et al. 2009a; and 
HD 209458b, Swain et al. 2009b). Both, CoRoT (Alonso et al. 
2009a, 2009b; Snellen et al. 2009) and Kepler (Borucki et al. 
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2009) are now also measuring occultations from space in the op- 
tical. 

In this paper, we present the first ground-based detection of 
//-band thermal emission from an exoplanet. Our target, WASP- 
19b (Hebb et al. 2010), is a Jupiter-mass planet in a 19-h orbit 
around a G8 V-type star, and is the shortest-period exoplanet cur- 
rently known. 



2. Observations 

We observed an occultation of WASP-19b by its // mag = 10.6 
host star with the cryogenic, near-IR imager HAWK-I (Pirard 
et al. 2004, Casali et al. 2006), mounted on Yepun of the Very 
Large Telescope (Paranal, Chile). HAWK-I is composed of four 
Hawaii-2RG chips, each measuring 2048x2048 pixels. Its pixel 
scale is 0.106"/pixel, providing a total field of view of 7'5 x 
7'5. Observations were obtained on 2009 May 03 from 23h06 
to 05hl7 UT. The broadband //-band filter was used {A e s = 
1 .620 //m, FWHM = 0.289 fim). The airmass ranged from 1 .09 to 
2.62 during the run and the transparency conditions were good. 
A total of 405 exposures, each comprising 10 integrations of 
1.26 s, were obtained. We used a pattern of six offsets, with the 
aim of producing an accurate sky map for each image from the 
neighbouring images. For a given offset, we kept the stars on 
the same pixels to minimise the effect of the small-scale spatial 
variations in the sensitivity of the detector chips. As the pointing 
was changed once early in the run, each star sampled a total of 
12 detector positions. To avoid saturation of the target and ref- 
erence stars, the telescope was heavily defocused, resulting in 
asymmetric stellar images with FWHMs of 13-28". 

In our analysis, we used only the images obtained with the 
Ql chip, which contained WASP-19 and several reference stars. 
After a standard pre-reduction (dark subtraction and flat-field di- 
vision), a localised smoothing was applied to each image: the 
count level of each pixel was compared to the median count 
level of the neighbouring pixels and, if the difference exceeded a 
threshold of 4 <x (or 30 cr for pixels belonging to stellar images), 
then the pixel's count level was set to the median of its neigh- 
bours. At this stage, a sky map was constructed for and removed 
from each image using a median-filtered set of the adjacent im- 
ages taken at different offsets. For each of the 12 offsets, aperture 
photometry was then performed using the IRAF/DAOPHOlQ soft- 
ware (Stetson, 1987). An aperture radius of 30 pixels was used. 
A similar flux extraction was also performed on the non-sky- 
subtracted images and found to provide a more reliable result. 
We attribute this to the incomplete removal of the large, asym- 
metric stellar images from the sky maps. We thus decided to use 
the fluxes extracted from the non-sky-subtracted images in our 
analysis, though we did measure the sky level in an annulus and 
subtract it from the stellar aperture. 

After a careful choice of reference stars, differential photom- 
etry was performed and a light curve was produced for each of 
the 12 offsets (Fig. 1). The light curve of one offset is signif- 
cantly poorer than the 1 1 others, probably because of a detec- 
tor defect in the image of either WASP-19 or a reference star, 
and we discarded this light curve from our analysis. For each of 
the 1 1 remaining light curves, the scatter is much larger during 
one half-hour period. This period corresponds to the minimum 
of the FWHM of the stellar images and to the maximum pixel 
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1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Top: The 405 HAWK-I measurements from the twelve 
pointing offsets. The light curve from each offset was given the 
same normalisation factor. Bottom: The 276 measurements re- 
maining after the rejection of one offset and an half-hour period 
that exhibited a non-linearity issue. The light curves were nor- 
malised on a per-offset basis. 



value (above 30kADU) for the target star. The larger scatter in 
this portion of the light curves is therefore probably caused by 
a non-linearity effect, and so we discarded these data. The 276 
measurements remaining (from an original 405 measurements) 
after rejection are shown in Fig. 1. 

Each light curve varies differently with time, and we found 
that these variations are strongly correlated with the background 
amplitudes. We chose not to simply detrend the light curves 
for external parameters and analyse the resulting corrected light 
curves. To avoid underestimating the error bars of our final pa- 
rameters, we included trend models in our global analysis (Sect. 



0. 



3. Data analysis 

To place as many observational constraints as possible on the 
occultation parameters, we performed a global analysis of our 
HAWK-I occultation photometry combined with the 34 radial 
velocities (RVs) and the FTS z-band transit light curve presented 
in Hebb et al. (2010). In addition, the SuperWASP transit epoch 
reported by Hebb et al. (2010) was used to constrain the orbital 
period of the planet. These data were adopted as input of the 
adaptative Markov-Chain Monte Carlo (MCMC) algorithm pre- 
sented in Gillon et al. (2009a, 2009b). This MCMC implementa- 
tion uses the Metropolis-Hasting algorithm (e.g., Carlin & Louis 
2008) to sample the posterior probability distribution of adjusted 
parameters for a given model. Our model was based on an oc- 
culting star and a transiting planet on a Keplerian orbit about 
their common centre of mass. We used a classical Keplerian 
model for the RVs obtained outside of transit (we discarded the 
single RV obtained during transit). To model the eclipse photom- 
etry, we used the photometric eclipse model of Mandel & Agol 
(2002), multiplied by a systematic effect model. 
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For the FTS transit, a quadratic limb-darkening law was as- 
sumed. Quadratic coefficients u\ and u 2 were interpolated from 
the tables of Claret (2004) for the Sloan z'-filter and for T eS 
= 5500 + 100 K, logg = 4.5 ± 0.2, and [Fe/H] = 0.02 ± 0.09 
(Hebb et al. 2010). We obtained u x = 0.250 + 0.015 and u 2 = 
0.304 + 0.008. We allowed the quadratic coefficients u\ and u 2 
to float in our MCMC analysis, using as jump parameters^ the 
combinations c\ — 2 x u\ + u 2 and c 2 - u\ - 2 x u 2 to minimise 
the correlation of the uncertainties (Holman et al. 2006). To ob- 
tain a limb-darkening solution consistent with theory, we added 
to our merit function the Bayesian penalty 

where c' j is the initial value deduced for the coefficient c, and cr c < 
is its error. 

For each light curve, the eclipse model was multiplied by 
a trend model to take into account known low-frequency noise 
sources (instrumental and stellar). For the FTS transit light 
curve, we modelled this possible trend as a quadratic time- 
dependent polynomial. Given the small number of points in the 
HAWK-I light curves obtained prior to the re-pointing of the 
telescope, we used a simple linear time-function as a trend model 
for them. For the other five HAWK-I light curves, our adopted 
trend model was a combination of a quadratic function of time 
and a quadratic function of the background amplitude. As ex- 
plained in Gillon et al. (2009b), the coefficients of the trend mod- 
els are not jump parameters of the MCMC, but are rather deter- 
mined by least squares minimisation at each step of the MCMC. 
We followed the procedure described in Winn et al. (2008) to 
check for correlated noise in the photometric time-series and to 
scale the error bars accordingly. For the RVs, the quadratic addi- 
tion of a jitter noise of 13 m s _1 was needed to obtain a residual 
rms in good agreement with the mean error in the measurements. 

The jump parameters in our MCMC simulation were: the 
planet/star area ratio (R p /R s ) 2 ; the //-band occultation depth; the 
transit width (from first to last contact) T^; the modified impact 
parameter b' — a cos i/R* (Gillon et al. 2009a); the orbital period 
P; the time of mid-transit TV, the two Lagrangian parameters 
e cos u> and e sin to, where e is the orbital eccentricity and a> is 
the argument of periastron; and K 2 = - e 2 /P 1 ^ (Gillon 

et al. 2009a), where K is the semi-amplitude of the radial stellar 
reflex velocity. We assumed a uniform prior distribution for all of 
these jump parameters. The merit function used in our analysis 
was the sum of both the^ 2 of each time-series and the Bayesian 
penalty presented in Eq. 1 . 

At each step of the MCMC, stellar mass was determined us- 
ing the stellar mass calibration relation of Torres et al. (2009), 
which was derived from detailed studies of detached stellar bi- 
naries. The parameters in this relation are log g, T e ff and [Fe/H], 
though, following the method of Enoch et al. (2010), we recast 
it in terms of stellar density, p*, instead of logg. Stellar density 
is constrained by observation, as it depends on the shape of the 
transit light curve and the eccentricity of the orbit, which is con- 
strained by the RVs and occultation timing. The values of T e ff 
and [Fe/H] were drawn from the distributions A^(5500, 100 2 )K 
and N(Q.Q2, 0.09 2 ), respectively, where N(fi, cr 2 ) is a normal dis- 
tribution with mean fi and variance cr 2 . To account for the uncer- 
tainty in the parameters of the stellar calibration law, the values 
of these parameters were randomly drawn at each step of the 

2 Jump parameters are the model parameters that are randomly per- 
turbed at each step of the MCMC. 



HAWK-I ff-band FTS z-band 




-10 -5 

Time from mid-transit [hours] 

Fig. 2. Top: Photometry of the planetary transit (FTS data, from 
Hebb et al. 2010) and the planetary occultation (HAWK-I data), 
folded on the transit ephemeris. The best-fit light curve model is 
superimposed. Bottom: Radial velocities (from Hebb et al. 2010) 
folded on the best-fit transit ephemeris. The best-fit RV model is 
superimposed. 



MCMC from the normal distributions presented in Torres et al. 
(2009). 

We first perfomed a single MCMC to assess the levels of cor- 
related noise in the photometry and jitter noise in the RVs, and to 
scale the measurement error bars accordingly. This MCMC con- 
sisted of 10 5 steps, with the first 20 % considered as its 'burn-in' 
phase and discarded. Five new MCMCs (10 5 steps each, with 
20 % burn-in) were then performed using the scaled error bars. 
The good convergence and mixing of these five MCMCs was 
verified using the Gelman and Rubin (1992) statistic (R parame- 
ter < 1 %). Finally, the inferred value and error bars of each pa- 
rameter were obtained from its marginalised posterior distribu- 
tion. Table 1 shows the deduced values of the jump and derived 
parameters. The best-fit planet model is shown in Figs. 2 and 3. 
As can be seen in Fig. 3 and Table 1, the occultation is clearly 
detected, with a depth of 0.259+° Jj^ %. The rms of the HAWK-I 
residuals is 0.12 %, and 0.047 % after binning in intervals of 10 
measurements. 

To investigate the robustness of the marginal detection of a 
non-zero eccentricity (e = 0.016+°^), we also perfomed an 
MCMC with eccentricity fixed to zero. To compare the eccentric 
and circular models, we computed the Bayes factor (e.g., Carlin 
& Louis 2008), which is the ratio of the marginal likelihoods 
(e.g., Chib & Jeliaskov 2001) of the two models. The deduced 
Bayes factor of 88 indicates that the data are more consistent 
with the eccentric model, though not decisively so. 



4. Discussion 

We have measured a decrease in the //-band flux from the 
WASP-19 system of 0.259+J}^ %, which we attribute to the oc- 
cultation of day-side planetary thermal emission by the host star. 
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Parameter (unit) 



Value 



Fitted model parameters (jump parameters ) 



b' = a cos i/R, 
T H (days) 
T (HJD) 
P (days) 

K 2 (ms- 1 days 1 ' 3 ) 
ecos a) 
esinw 

//-band F p /F, 



01982+ 000050 

-0.00049 

631+ 0026 

06506+ 000072 
u.uujou 00007 , 

2454776.91490 ±0.00019 

0.7888393 ± 0.0000045 

236.9!^ 

0069+ 00024 

u.uuu5> 0027 

0.005 ±0.021 

00259+ 00046 

—0.00044 

0.810 ± 0.030 
-0.357 ± 0.022 



Dependant parameters deduced from the above 



K(msr l ) 

h 

"transit 
^occultation 
^occultation (HJD) 

a(AU) 

in 

e 

con 

M, (M ) 
R, (Re) 
P, (Po) 
ill 
«2 

M„ (Mj) 
R P (Rj) 
P„ (Pj) 



256.5 ± 6.0 

°- 628 -a029 

633 +OOM 
_ 0034 

2454777. 31 28 1!™™ 3 
01659+ 000040 

vJ. l —0.00039 

80.07!"? 2 
016+ 0015 

-0.007 
a a +44 

J -101 

fl Q77+0.073 
"• y ' ' -0.067 

975 +a038 

7 -0.037 

1 054 +0095 
253 +a015 

0.3049 ± 0.0083 



1.166! 
1.338! 
0.486+ 



Table 1. WASP-19 system parameters and l-cr error limits. 
1 Constrained parameter; see text for details. 



Our global solution suggests that the planet's orbit may be 
non-circular, the value of e cos to being non-zero at the 2.6-cr 
level; the constraint on e sin to is far weaker. The measurement 
of additional occultations and RVs will constrain e cos to and 
e sin to more strongly. A non-circular orbit may indicate that 
tidal heating is responsible for WASP-19b's inflated radius (e.g., 
Ibgui et al. 2009). 

The measured occultation depth corresponds to an //-band 
brightness temperature of 7# = 2580 + 125 K. To calculate Th, 
we defined the product of the planet/star area ratio and the ra- 
tio of the bandpass-integrated planetary to stellar surface photon 
fluxes, corrected for transmission^], to be equal to the measured 
occultation depth (e.g., Charbonneau et al. 2005). We assumed 
the planet to emit as a black body and, for the star, we used a 
model spectrum of a G5V star (Pickles 1998), normalised to re- 
produce the integrated flux of a black body with T e s = 5500 K 
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4955.55 4955 6 4955.65 4955.7 

BJD - 2454000 

Fig. 3. HAWK-I photometry corrected by the best-fit trend 
model for each time-series. The best-fit occultation model is su- 
perimposed. Top: unbinned. Bottom: binned every 10 measure- 
ments. Our global solution indicates a mid-occultation time of 
BJD=2 454 955.5905. If a circular orbit is assumed, the transit 
ephemeris of Table 1 predicts a time of BJD=2 454 955.5870, 
earlier by 0.0035 d (=5 min =0.35x a minor tic increment). 



(Hebb et al. 2010). The uncertainty in Th only takes into account 
the uncertainty in the measured occultation depth. 

In Fig. 4, the measured //-band flux density ratio is com- 
pared to a model atmosphere spectrum of the planet (Barman et 
al. 2005), which was based on parameter values from Table 1. 
The model is cloud-free and adopts the nominal day-side sce- 
nario, in which the incident stellar energy remains entirely on 
the day-side. Given the planet's small orbital distance from the 
G8V-type host star, the irradiation is quite strong (r eq « 2400 K, 
day-side only) and, in these models, it is common to see a nearly 
isothermal photosphere with a large temperature inversion at 
low pressures. In the model explored here, the inversion layers 
have temperatures exceeding 3000 K. Despite these high tem- 
peratures, the inversion does not extend deep enough into the 
photosphere to match the observed //-band flux; the model un- 
der predicts the flux by roughly a factor of two. 

To reconcile our model's prediction with the occultation 
measurement, one would need to increase the luminosity of the 
star by a factor of two, for example, by increasing the stellar tem- 
perature by 1000K or increasing the host star radius by 0.4 R G 
(or some combination of these). These values of the stellar pa- 
rameters are excluded by observations (Hebb et al. 2010; this 
work). The most probable explanation of the poor match to the 
planet flux is therefore related to the model atmosphere. 

It is possible that the day-side atmosphere is much closer to 
pure radiative-convective equilibrium and experiences little day- 
side redistribution. This zero-redistribution scenario has been 
shown to result in a hotter substellar point and deeper occul- 
tations than the uniform day-side-only model (Barman et al. 
2005). As more data are obtained, in particular at other near- 
IR wavelengths, we will explore a wide variety of atmospheric 
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Fig. 4. Synthetic planet/star flux ratios of a cloud-free model that 
assumes solar abundances and no redistribution of incident stel- 
lar flux to the night-side (solid line). The effective temperature 
of the star was taken to be T eS = 5500K (Hebb et al. 2010); 
all other parameters are taken from Table 1 . The flux ratio for a 
black-body planet (at 2400K) is indicated by a dotted line. Solid 
symbols are the //-band integrated flux ratios; symbols without 
error bars correspond to the model points. The HAWK-I //-band 
filter response curve is also shown. 



phenomena (e.g., clouds, depth-dependent energy redistribution, 
and photochemistry). 

Swain et al. (2010) found a strong, unexpected 3.25-fj.m 
emission feature in the day-side spectrum of HD 189733b, which 
they attributed to non-local-thermodynamic -equilibrium (non- 
LTE) CH4 emission. Future work should explore whether non- 
LTE chemistry is responsible for the unexpectedly high //-band 
flux of WASP- 19b. For now, our measured //-band flux presents 
an interesting puzzle for irradiated atmospheric models. 

Acknowledgements. M. Gillon acknowledges support from the Belgian Science 
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